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Overlap with the separable state is introduced in this paper for the purpose of characterizing the 
overall correlation in many-body systems. This definition has clear geometric and physical mean- 
ing, and moreover can be considered as the generalization of the concept- Anderson Orthogonality 
Catastrophe. As an exemplification, it is used to mark the phase transition in the Dicke model 
for zero and finite temperature. And our discussion shows that it can faithfully reflect the phase 
transition properties of this model whether for zero or finite temperature. Furthermore the overlap 
for ground state also indicates the appearance of multipartite entanglement in Dicke model. 
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I. INTRODUCTION 

Correlation in condensed matter systems predominates 
the understanding of many-body effects. Fundamentally 
one can define different correlation functions for describ- 
ing the unusual connections in many-body systems. For 
instance it is general to introduce the order parameter 
for the description of phase transitions induced by lo- 
cal perturbation, and furthermore to classify the diverse 
phase transitions by scaling the singularity of correlation 
functions with the universal critical exponents. That is 
so called Landau-Ginzburg-Wilson(LGW) paradigmjf]. 
However the situation becomes different for the strongly 
correlated electronic systems. The quantum Hall ef- 
fect appearing in two-dimensional electron gas with high 
magnetic field shows the distinct features not captured by 
LGW paradigm. Instead topological order consequently 
is defined to describe the underlying symmetry in quan- 
tum Hall systems, which is distinct from the notion of 
spontaneously broken symmetry 0. 

Recently the extensive researches of quantum entan- 
glement in condensed mattered systems show the po- 
tentiality that quantum entanglement would act the 
universal description for many-body effects^, E[. Es- 
pecially some general conclusions have been obtained 
about the connection between quantum entanglement 
and quantum phase transition in many-body systems. 
The concurrence, a measurement of two-party entangle- 
ment, has been shown to behave singularly at the crit- 
ical point of one-dimensional spin-half XY model, and 
the critical exponents can also be obtained by scaling 
this singularity [(| Q. Furthermore the block entangle- 
ment entropy has been shown to display the logarithmi- 
cal divergency with the block size at the critical points, 
and the scaling factor is directly related to the central 
charge of the conformal field theory Q. Moreover the 
universal area law for the entanglement entropy has also 
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been constructed exactly in one dimensional spin-chain 
systems Q, and the similar behavior for single-copy en- 
tanglement is also founded^. Recently the entangle- 
ment spectrum has been defined to obtain the general 
information about the many-body systems @-|H|- As for 
quantum Hall systems, it is shown that the scaling be- 
havior of entanglement entropy is directly related to the 
quantum number, which is used to characterize the topo- 
logical order [l2j • And entanglement spectrum can also 
be used to detect the non-local features of quantum Hall 
systems 

Although these important progresses have been made, 
there are a few exceptions that lead to the suspicion of the 
validity of quantum entanglement as an universal descrip- 
tion for many-body effects. Entanglement entropy some- 
times provides ambiguous information about the phase 
transitions in higher dimensional many-body systems 
Even for one dimensional systems, it cannot present the 
complete information in some situations. As an exam- 
ple, the recent studies show that the block entanglement 
entropy for the Valence-Bond-Solid(VBS) state of inte- 
ger spin seems unsensible to the degeneracy manifested 
by the underlying topological symmetry and also does 
not display the dependence on the parity of spin num- 
ber s, which however both can be manifested clearly by 
introducing string order parameter [13|. As for quantum 
Hall systems, the entanglement entropy and entangle- 
ment spectrum have also been shown the limited ability 
of identifying the topological orders [Tl|. 

In my point, this defect would attribute to the trace- 
out of the superfluous degrees of freedom when one ob- 
tains the reduced density matrix. And some information 
for the global features in many-body systems is inevitably 
lost. This point has been exemplified in a recent paper of 
our group[14j], in which the geometric entanglement(GE) 
as a measurement of multipartite entanglement is calcu- 
lated for VBS state. The interesting result in this paper 
is that GE displays two different scaling behaviors de- 
pendent on the parity of spin number s, and the global 
GE is divergent linearly with the particle number. 

Through this short introduction, it seems promising to 
measure multipartite entanglement in order to obtain the 
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complete information for many-body effects. Recently 
some efforts have been made in this direction. The con- 
nection of multipartite entanglement and quantum phase 
transition has been discussed in some special models [T^- 
lvi\ . However the crucial obstacle for further development 
is the absence of the unified understanding of the multi- 
partite entanglement [l8l [l9j . Whereas the maximally en- 
tangled state can be defined unambiguously for bipartite 
systems, what is the maximally entangled state for mul- 
tipartite systems is unclear until now[l9j. Fortunately 
it is well accepted that the fully separable state can be 
defined as 
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where TV is the particle number, and pi denotes the 
common probability with which the single-particle state 
Pn\n = 1,2,- •• ,N) happens. With respect of this 
point, geometric entanglement(GE) is introduced first by 
Shimony for pure bipartite state[20| and generalized to 
the multipartite case by Carteret et aZ.[2l[, Barnum and 
Linden [22T |. Wei and Goldbart [23j . and to the mixed 
state by Cao and Wang [2J|. GE is a genuine multipar- 
tite entanglement measurement. The main idea of GE 
is to minimize the distance D between the state |*) to 
be measured and the fully separable state |$) in Hilbert 
space, 



D = minllH*) - l$)H 2 }- 
{|t)) 



(2) 



For the normalized I*) and the evaluation of D is 
reduced to find the maximal overlap 23] 



A(*) = max|($|*)|. 
{I*}} 



(3) 



Geometrically Ad*)) depicts the overlap angle between 
the vectors |*) and |$) in Hilbert space. Then the larger 
Ad*)) is, the shorter is the distance and the less en- 
tangled is |*). But the optimum is in general a forbid- 
den task, not spoken for mixed state, and the analyt- 
ical results can be obtained only for some very special 
cases [H, [13, [HI- Recently many efforts are devoted to 
the reduction of the optimum and some interesting re- 
sults are obtained pBj. 

Given this difficulty, we introduce another different 
quantity in this paper to capture the overall correlation in 
condensed matter systems, i.e. the overlap with a special 
fully separable state. The starting point is still to find the 
minimal distance between the state to be measured and 
a special fully separable state defined in the next section. 
In contrast to GE the optimum can be reduced by utiliz- 
ing the geometric property of the overlap, and this over- 
lap has very clear physical meaning, whether for pure or 
mixed state. In Sec. II, the definition of this overlap is in- 
troduced, and the differences with several known similar 
definitions are clarified. Furthermore we point out that 
our definition is connected intimately with the concept of 



Anderson Orthogonality CatastrophefAOC)[26l [27j. As 
an illustration of the validity of our definition, the col- 
lective phase transition appeared in Dicke model is dis- 
cussed by this quantity in Sec. III. Multipartite entangle- 
ment in this model is also studied for displaying the po- 
tential connection between this overlap and multipartite 
entanglement. Finally conclusions and further discussion 
are presented in Sec. IV. 



II. OVERLAP WITH FULLY SEPARABLE 
STATE 

Similar to the introduction of GE, our starting point 
is also to find the minimal distance D between the fully 
separable state p sep and the state p to be measured. 
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Generally this minimal distance is still decided mainly by 
the maximal overlap 



A = max Tr[pp sep ] 



(5) 



The density matrix can also be written as the Bloch 
form 



p = (i+J2 r * A *)M 



(6) 



where d denotes the dimension, Ai is the generator of 
SU(d) group and {ri} is so called Bloch vector [30]. Thus 
A has clear geometric meaning which depicts the mini- 
mal overlap angle 6 between the Bloch vectors {r^} and 
{fi}sep in the Bloch- vector space, i.e., 



m & xTr{pp sep ] = -(l + \{r l }\\{r i } 
{p°'p} d 
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Two limit cases are beneficial to the understanding of the 
physical meaning of 8. For cos 9 = 1, the overlap is max- 
imal and p and p sep share the same physical characters 
since Bloch vector {r{} is the reflection of the intrinsic 
symmetry in the systems [3(| . While for cos 6 = — 1 one 
has minimal overlap, and p sep shows distinct properties 
from that of p. 

In contrast to the Bures fidelity |3lj|. the overlap A have 
clear geometric meaning whether for pure or mixed state, 
as shown above. Furthermore by this geometric meaning, 
the optimal procession can be reduced to find the fully 
separable state p sep sharing the same physical properties 
with p (see Appendix A for a proof). Moreover this def- 
inition is more popular than Eq.@. First A comes back 
to the form Eq.© for pure states. Second Eq.© includes 
the case when one state is pure and the other is mixed. 
This situation always happens as exemplified in Ref. [24j , 
but is not covered in the original discussion [2 3]. Thirdly 
for mixed state the geometric characters of GE becomes 
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ambiguous because of the convex roof construction [23j, 
while the geometric meaning of A is clear whenever for 
pure or mixed state. 

With these advantages, the evaluation of A however is 
difficult for mixed state p sep since there are infinite possi- 
bilities of the decomposition for p sep . Recently we note a 
popular concept in condensed mat ter phys ics- Anderson's 
Orthogonality Catastrophef AOC) |26l |27I |. which refers 
to the vanishing of the overlap between the many-body 
ground states with and without the potential as a power 
law in the number of particles in the systems. AOC is 
defined as 



(8) 



where \& p ) and |<&) correspond respectively to the many- 
body states with the potential and the state described en- 
tirely in terms of free plane waves, including the ground 
state of the unperturbed system 26]. Anderson proved 
that the overlap A approached to be zero under thermo- 
dynamic limit N — > oo even for very weak potential, that 
means that |$) is orthogonal to |$ p )|2(| and the transi- 
tion between the two different states is forbidden. It is 
the physical meaning of catastrophe [27l| ■ As claimed in 
Ref. [26[ , it becomes impossible because of the appearance 
of catastrophe to find the characters for many-body sys- 
tems by adiabatically imposing the potential and observ- 
ing the response, since the significant changes in many- 
body systems can be induced even for infinitesimal per- 
turbation. AOC presents an understanding of a number 
of Fermi-edge singularities, e.g., in the Kondo effect (28l| 
or in the X-ray edge problem[29], for which the local 
singularity has an overall effect on the property of the 
many-body systems. With these points AOC manifestly 
shows that the correlation in many-body systems can be 
constructed simultaneously whenever the interaction ap- 
pears, and thus can be used to give a full description 
of correlation in many-body systems. Furthermore the 
important feature is that this prohibition can be con- 
quered by the symmetry-breaking process as exemplified 
by the observation of the X-ray absorption in the electron 
gas [27J, which means that AOC can also be used to char- 
acterize the phase transitions induced by the symmetry 
breaking process. In a word AOC presents an compre- 
hensive description of the many-body effects. 

This crucial observation enforces us to define the fol- 
lowing fully separable state for N parties, 
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which represents the many-body state without potential, 
compared to the state |$) in Eq.©. And then we can 
define the overlap with fully separable state p s to capture 
the overall correlation in many-body systems, 



A = maxTr[pp s 
{p 3 } 



(10) 



1. A has clear geometric meaning, which depicts the 
minimal overlap angle between the Bloch vectors 
{ri} and {ri} s . And for pure states, it returns to 
the original definition Eq.Q of GE. 

2. By this geometric meaning, the optimal process in 
A can be reduced to find the fully separable state 
p s sharing the same physical features with p. 

3. A can be regarded as the generalization of AOC to 
mixed state, and can faithfully reflect the overall 
correlation in many-body systems. 

It should emphasize that this definition does not try to 
present a complete measurement of the multipartite en- 
tanglement. Instead our purpose is to find an universal 
method to characterize the overall correlation in many- 
body systems, whether quantum or classical. However 
this definition is also meaningful to find the unified un- 
derstanding of multipartite entanglement in many-body 
systems. As shown in the next section, A indeed presents 
the interesting information for the phase transition in 
Dicke model. And moreover the connection between A 
and multipartite entanglement in Dicke model has also 
been discussed in Sec. III. Additionally in contrast to the 
recent interest in the fidelity for many-body systems |32j|. 
A does not serve for the state discrimination. 



III. EXEMPLIFICATION: PHASE TRANSITION 
IN THE DICKE MODEL 

In order to demonstrate the generality of this defini- 
tion, the phase transition in Dicke model is discussed by 
A in this section. Dicke model describes the dynamics 
of N independently identical two-level atoms coupling 
to the same quantized electromagnetic field (33[. Due to 
the presence of dipole-dipole force between atoms^ Dicke 
model shows the normal-superradiant transition[34l|. 

Dicke model is related to many fundamental is- 
sues in quantum optics, quantum mechanics and con- 
densed matter phy sics, such as the coherent sponta- 
neous radiation 3J| , the dissipation of quantum system 
[35j . quantum chaos[36j and atomic self-organization in 
a cavity [13] ■ The normal-superradiant transition have 
been first observed with Rydberg atoms [3{|, and re- 
cently in a superfluid gas coupled to an optical cavity [37j 
and nuclear spin ensemble surrounding a single photon 
emitter [381] . Quantum entanglement in Dicke model has 
also been discussed extensively in[4(| l4l| . Furthermore 
Dicke model is also related to the issues of how the 
opened multipartite systems is affected by the environ- 
ment and the robustness of multipartite entanglement [42l| 

The Hamiltonian for single-model Dicke model reads 
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This definition is the main contribution in this paper, and 
has some distinct advantages, summarized as following 
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where J z = J2iLi °f /2 an d J± — Ya=i a f are tne co ^" 
lective angular momentum operators. At zero tempera- 
ture, the normal-superradiant transition happens when 
A = A c = ^/lJuJq/2. For finite temperature, the critical 
temperature is decided by the relation[43] 



ujq tanh(/3 c w/2) 
2A 1 tanh(/3 c o; /2) ' 



(12) 



An intrinsic property of Dicke model is the parity sym- 
metry, 



[H,U] = 0, 

Yl = giw(a^ a+J z + - 



(13) 



Moreover Eq. (|lip is obviously permutation invariant by 
exchanging any two atoms. 

With these information, the overlap A for the Dicke 
model is studied explicitly for zero and finite tempera- 
ture in the following two subsections. Some interesting 
features of A are displayed. 



A. Zero Temperature 



With respect of the permutation invariance of atoms 
in Eq. (flTT) , it is convenient to introduce the Holstein- 
Primakoff(HP) transformation 
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with bosonic operator b^\ Semiclassically there is a 
ground state with J z = —N/2 for Dicke model under 
thermodynamic limit N — > oo. Hence it is reasonable to 
adopt the low-excitation approximation at zero temper- 
ature, and then one obtains two effective Hamiltonians 
for different regions of A (refer to [44| for details) 



J 



= wa t a + wo& t & + A(a f +a)(6 f +6) - yw , A < A c ; 
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(15) 



if W and can be diagonalized readily by transform- 
ing into phase space, and then one has the diagonalized 
form El 



(16) 



where the forms of w 1 ( 2 ) and C]Y2) are dependent on A > 
A c or A < A c [44j]. Then the ground state can be written 
as 



Is) = |s)s 



(17) 



where 



'1(2) 



denotes the vacuum state for mode uju^)- 



Furthermore the average spin along z direction per atom 
shows distinct values across the phase transition point, 



TV 



■s, A < A c ; 

A 2 

2A^' A > A c , 



(18) 



which then can act the order parameter. It is obvious 
that a macroscopic number of atoms are excited for A > 
A c , which is so called superradiant phase, while for A < A c 
it is normal phase. 

With these information, we are ready to evaluate A. 
Our focus is mainly on the atom system. Then the crucial 



step is to decide the fully separable state p s for atom 
system. As mentioned in Sec. II and proved in Appendix 
A, the optimum process in Eq.(|4|) can be reduced to find 
p s sharing the same global features with the ground state 
Eq. (|17p . First with the requirement of the permutation 
invariance of atoms in Dicke model, the single atomic 
state should have the same form in p s , i.e. pi — g,i — 
1, 2, • • • ,N, and then 



P = Q 



(19) 



Second the parity symmetry for Dicke model is reduced 
for single atom state p as 



[e" J %p s }=0 
=>[e*"",Q]=0. 

Thus one has under a z representation 
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FIG. 1: The overlap A with fully separable state p s vs. the 
coupling A at zero temperature, uo = oj = 1 has been chosen 
for this plotting, and the critical point is A c = 0.5 in this case. 
The inset is a plotting for the purity of the reduced density 
of atomic freedom under N — > oo. 



Finally with requirement of Eq. (|18 
p s can be uniquely determined as 



1/2- 



■U*L. Thus 




(22) 



It should point out that the procedure for the determi- 
nation of p s is popular whether for zero or finite temper- 
ature. 

For the evaluation of the overlap A, it should note that 
p s can be rewritten as the following contract form under 
the J z representation, 
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where atC/c denotes the binomial function, and |n 
presents the state for which n particles is spin-up and the 
other is spin-down. Together with the HP transformation 
Eq. ([T4| . it is obvious 
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Then A can be evaluated easily under this representation. 

In FigQ] the overlap A with p s is plotted. At nor- 
mal phase (A < A c ), one has (J z )/N — —0.5, a — 0, and 
then p s is the fully separable pure state. A is determined 
mainly by the first diagonal element of the reduced den- 
sity matrix of atom system in this special case. While 
for A > A c , A shows a sudden rising and then decreases 
with the increment of A, and tends to be steady with 
A — > oo. Moreover under N — > oo, A tends to be vanish- 
ing. Then two different phases can be clearly identified 
by evaluating A. 



Some intricate features of the phase transition can be 
disclosed by A. For normal phase A < A c , it is known 
that the atom system becomes entangled with the elec- 
tromagnetic field, and attains the maximal value at the 
critical point [40j. The entanglement leads the state of 
atom system to be mixed, and the purity of its reduced 
density is decreased as shown by the inset in Fig[TJ At 
the same time the pairwise entanglement between two 
any atoms is also raised mediated by their couplings to 
the electromagnetic field, and the atoms become corre- 
lated with each other [4(| • These intrinsic properties can 
be captured by A at the same time. For normal phase 
p s is pure and fully separable. Thus the decrement of A 
reflects the fact that the atoms become correlated with 
each other, and attains the maximal correlation at the 
critical point, at which A has minimal value. Further- 
more since there is no interaction among atoms, the only 
reason for the construction of correlation in atoms is the 
couplings to the same electromagnetic field, which just 
induces the state for atom system to be mixed. This 
feature can also manifested by the decrement of A with 
respect that p s is pure. 

For superradiant phase A > A c , it is known that the en- 
tanglement between the atoms and electromagnetic field 
decreases monotonously to a steady value with the incre- 
ment of A, while the pairwise entanglement in atoms dis- 
appears asymptotically [40| . Contrastably the purity for 
the state of atom system has a sudden increasing closed 
to A c and then decreases to a steady value, as shown by 
the inset of Fig[TJ The two different behaviors can also 
be captured by A. Similar to the behavior of the purity 
of the state for atom system, A has also a sudden arising 
closed to A c and then decreases to a steady value with 
the increment of A. With respect that p s is mixed in this 
case and its purity is monotonically decreased with the 
increment of A, the abrupt increment of A means that 
the sudden recovery of the purity of atomic system is at 
the expense of the reduction of correlation in atoms. It 
is obvious from FigQ] that A tends to be zero with the 
increment of N for large A. However the vanishing of A 
cannot attribute to the mixedness of p s since the steady 
value of A for finite N is always bigger than the maximal 
mixedness 1/N, manifested by Fig[TJ This feature means 
that the correlation in atoms still exists. Since the pair- 
wise entanglement of atoms is known to be vanished in 
this limit [401] . the correlation in atoms must be global. 

The scaling behavior of A near the critical point show 
some interesting features. At the normal phase (A < A c ), 
one has for uj = luq = 1 



A = 



2 3/2 (1 _ 4A2) l/4 



[1 + 3Vl - 4A 2 + 0.5(Vl + 2A + Vl - 2A) 3 ] 



(25) 



Similar to the method in Ref.|17j. one can define the 
globe overlap — In A to measure the atomic correlation 
in Dicke model. It is obvious that the globe overlap is 
mainly determined by (1 — 4A 2 ) 1 / 4 near A c = 1/2, and 
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then 



-lnA~-iln(l-£) 



(26) 



which is same to the scaling behavior of multipar- 
tite entanglement in the Lipkin-Meshkov-Glick(LMG) 
model[I3]- This result is not strange since Dicke model 
and LMG model belong to the same universality class. 
However it strongly implies that A could be correlated 
directly to the multipartite entanglement in Dicke model. 
As shown in Sec. IV, the atom system indeed displays the 
multipartite entanglement in this case. 



B. Finite Temperature 

At finite temperature the phase transition is induced 
by thermal fluctuation. In order to determine the crit- 
ical temperature, the general method is to evaluate the 
partition function z. In Ref . (43| . z has been obtained 
analytically 
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and the critical temperature is determined by Eq. (ll2[) . 
For w — Do, it is reduced to T c = 2 A 2 /fegwo- With 
the same trick used in [43|, the overlap A can also be 
written analytically as (see Appendix B for the details of 
calculation) 
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As shown in Fig[2l A can clearly detect the phase tran- 
sition by its abrupt variance closed to the critical line. 
With respect that p s is mixed and fully separable in this 
case, A reflects that the correlation in atoms exists even 
for finite temperature. However this type of correlation is 
obviously induced by the thermal fluctuation, and thus 
is incoherent in contrast to that for zero temperature. 
This difference will become clear if one focuses on the 
multipartite entanglement of atoms in the next section. 



IV. MULTIPARTITE ENTANGLEMENT IN 
DICKE MODEL 

Another interesting aspect for Dicke model is the mul- 
tipartite entanglement in atoms. Since all atoms simulta- 
neously couple isotropically to the same electromagnetic 
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FIG. 2: Overlap A with fully separable state p s vs. the cou- 
pling A and temperature T. ujq = ui = fcs = 1 and TV = 100 
has been chosen for this plotting. The two 3-dimensional fig- 
ures are the same plotting with different plot-ranges for clar- 
ity. In the contour plotting the red dashed line corresponds 
to the critical line T c — 2\ 2 /kstoo, and because of the rapid 
decay of A only finite range of its values is shown for this 
contour plot. 



field, then it is expected that the multipartite correlation 
of atoms could be readily constructed in this case. 

However the measure of multipartite entanglement is 
a difficult task in general, especially for mixed state. An 
indirect way of resolving this difficulty is to find the char- 
acters uniquely belonging to the fully separable state 
Eq.([T]), and the violation of these properties implies the 
appearance of multipartite entanglement. Spin squeez- 
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FIG. 3: Eqs. (|30|) vs. the coupling A at zero temperature. 
u)q — cj — 1 and TV = 100 have been chosen for this plotting. 
The labels z — x — y and y — z — x denote the sequence and 
values of a — ft — 7 in corresponding inequalities. 

ing is one of the most successful approaches to the multi- 
partite entanglement in this way[45j. Recently G. Toth, 
et. al. provides a series of inequalities about spin squeez- 
ing to identify the multipartite entanglement in collective 
models [4^. 

{Jl) + {Jl) + {Jl) < (29a) 
A 2 J x + A 2 J y + A 2 J z > f, (29b) 
(J 2 a ) + {4)-f < (iV-l)A 2 J 7 , (29c) 
(7V-1)[A 2 J Q + A 2 J,] > (J 2 ) + A(A^) ; (29d) 

where ct,/5, 7 adopt the all permutation of x, y, z, and 
A 2 J a = (J 2 ) — (J Q ) 2 . The violation of any one of these 
inequalities implies the appearance of entanglement [46]. 
With respect of the limit large N, these inequalities can 
be rewritten as 

M( J ^ + ( J y) + ( J ")) < h (30a) 
^(A 2 J, + A 2 J, + A 2 J,)-^ > 0, (30b) 

h£ - M( J l) + (4)) + m > °> (30c) 

^(A 2 J a + A 2 J^)-^-i > 0, (30d) 

in which -^j(J 2 } and ^jA 2 J a are equivalent to evaluate 
the average {{J a /N) 2 ) and A 2 {J a /N) = ({J a /N) 2 ) - 
(J a /N) 2 . For large N, these inequalities have nontrivial 
result since the average magnetization per particle and its 
fluctuation still have nonvanishing values. It should point 
that Eq. (|30ap is obviously satisfied for arbitrary state. So 
the following discussion is mainly about Eqs. ([30"l) fb-d). 

A. Zero Temperature 

The evaluations of (J a /N) and ((J a /N) 2 ) 
can be implemented readily through Bogoliubov 
transformation Our calculations show that Eq. (|30bp 



is always satisfied at both normal and superradiant 
phases. In FigJ31 several situations for Eqs. (j3"0")) have 
been plotted with limit N —> 00, and the others can be 
proved to be bigger than zero. The violation implies 
that the atoms should be entangled. Moreover since the 
pairwisc entanglement between atoms is known to be 
vanished with increment of A [13], this entanglement is 
sure to be multipartite. Furthermore there is also a sud- 
den increment closed to the critical point, similar to the 
behavior of A shown in Fig[T] This feature means that 
there is a sudden reduction of the correlation of atoms, 
and A can also be used to detect the entanglement of 
atoms in Dicke model at zero temperature. 



B. Finite Temperature 

At finite temperature, the evaluations of (J a /N) and 
((J a /N) 2 ) can adopt the same trick used in Ref.[43|] (also 
shown in Appendix B). In FigsJU Eqs. ([30)) fb-d) have 
been plotted with all possible permutation of x, y, z. It 
is obvious that all inequalities are satisfied simultane- 
ously, and then one can conclude that there is no quan- 
tum entanglement of atoms in this case. This result is 
not surprising since the thermal fluctuation is dominant 
at finite temperature, and is considered to be incoherent. 
Although the absence of quantum correlation, the cor- 
relation induce by thermal fluctuation predominates, as 
shown in Figs 12] by A, which means that A can also be 
used to detect the thermal correlation. 



V. CONCLUSIONS AND FURTHER 
DISCUSSION 

In this paper, the overlap A with a special fully sepa- 
rable state defined in Eq. (fTU)) is introduced, in order to 
capture the overall correlation in many-body systems, 
whether quantum or classical. A has clear geometric 
and physical meaning shown in Sec. II. With these fea- 
tures, the optimum process in the definition of A can 
be reduced to find the fully separable state p s defined in 
Eq.®, which shares the same physical properties with 
the state to be measured. Importantly A can be con- 
sidered as the generalization of the concept Anderson's 
Orthogonality Catastrophe [l^, [27|, which is critical for 
the understanding of some effects in condensed matter 
physics. This important connection displays the popu- 
larity of A to detect the global correlation in many-body 
systems. And as an exemplification the phase transition 
in Dicke model has been discussed by A. 

As shown in Sec. Ill, A unambiguously depicts the 
phase transition features and the global correlation in 
Dicke model, whether for zero or finite temperature. 
At zero temperature, A displays the distinct behaviors 
across the critical point. Furthermore with the informa- 
tion of p s , A predicts the appearance of the multipartite 
entanglement in atom system, as verified in Sec.IVA. 
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Eq.(30b) 




FIG. 4: Egs. (|30[) vs. the coupling A and temperature T. 
luq — uj = ks = 1 and N = 100 has been chosen for these 
plots. The labels x-y-z,z ~ x — y and y — z — x denote the 
sequence and values of a — f3 — 7 in corresponding inequali- 
ties. The similarity among several plots is because their corre- 
sponding inequalities respectively would become closed under 
limit N — 7- 00. 



As for finite temperature, A can still be used to mark 
the phase transition in Dicke model, as shown in FigsH 
It displays the sudden variance at the critical line decided 
by the temperature T and the coupling A. An intricate 
feature appears when T — > 00. It is believed that all 
atoms would becomes independent in this case and can 
be considered as the ideal system[47]]. Quantum mechan- 
ically, the state in this case can be described by a fully 
separable state. Whereas A approaches zero as shown 
in Figs.@, and the non-zero A appears only at interme- 
diate temperature, shown in Figs. ([5]). This phenomenon 



implies that the correlation in atoms would exist even 
under high temperature. Moreover under the J z repre- 
sentation, the dimension is proportional to the atomic 
number N, and the value of the overlap shown in Figs. 
([2]) has exceeded greatly the limit by N. Thus this phe- 
nomenon cannot attribute to the mixedness of the state 
for atom system. Unfortunately we do not know how to 
understand these two different features. 

Although A cannot present a complete measurement 
of multipartite entanglement, it has been shown the in- 
~^timate connection to the quantum entanglement in some 
f special cases, such as the discussion for Dicke model at 
zero temperature in this paper. From this discussion A 
7' would presents a complete description for the global cor- 
relation in many-body systems, whether quantum or clas- 
sical. So it is not surprising that A can be used to identify 
the quantum entanglement in some special cases. How- 
ever it is difficult to answer the question what the general 
relation between A and quantum entanglement is since 
the absence of the unified understanding of multipartite 
entanglement. This point will be studied in the future 
publication. 

Appendix A: Find the nearest p s for a definite p 

For two arbitrary density matrixes p\ and p 2 , they al- 
ways has the following decompositions simultaneously 

Pi = ^P ( n\ n )n( n \ 

n 

Pi = XX^H^H, (31) 

m 

1(2) 

where denotes the probability that the system is 

I being in the state |n(m))i( 2 ). It should emphasize that it 
is unnecessary for the states labeled by different n or m 
to be orthogonal with each other. And thus the decom- 
positions above can always be realized at the same time. 
Then the overlap between p\ and P2 reads 

Tr[p lP2 ] = J2p { n ) P% ) \2{rn\n) 1 f. (32) 

m,n 

Obviously the maximization of overlap is dependent on 
the inner product |2(w|n)i| 2 . It is well known that for 
two different states \v) and \w) their inner product is 
bounded by Cachy-Schwartz(CS) inequality, i.e., 

\(v\w)\ 2 < (v\v)(w\w), (33) 

where the equality occurs if and only if the two vec- 
tors \v) and \w) in Hilbcrt state arc linearly related, i.e. 
\v) = c\w) for some scalar c. The important point for 
this condition is that c is not necessary a constant, for 
which \v) and \w) become physically identical, and CS 
inequality has trivial consequence. Thus 

Tr[p lP2 ] < J2plPpin ) i(n\n) 12 (m\m) 2l (34) 

m,n 
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where the equality occurs if and if only for arbitrary |n)i 
and 1 771)2 they are still linearly related. But in this case 
the scalar c has to be dependent on both n and m, i.e. 
c = c mn , which means that any (77)1 have to be linearly 
related to all 1 777)2- An interesting consequence for this 
condition is [pi,p2] = 0, which means that pi and P2 
share the same set of eigenvectors, and thus they share 
the same global symmetry and belong to the same space. 

This conclusion is not strange if one notes that the 
overlap between two matrixes is mainly determined by 
the inclusion relation of the spaces decided by the ma- 
trixes. As an example let consider two matrixes belong 
to two completely different spaces. And then the overlap 
must be zero since mathematically the intersection of the 
two spaces is null and there is no crossing items between 
the two matrixes. Comparably if one space is the sub- 
space or equivalent to the other space, the overlap then 
is nontrivial generally since the two matrixes belong to 
the same space. Hence in order to find the maximal over- 
lap between two matrixes, it is also necessary for the two 
matrixes to be in the same space. From physical point, it 
means that the two operator is necessary commutative. 
Furthermore it is easy to understand why the maximal 
GE for pure entangled state always happens for purely 
separable state. 

As for the determination of p s in Eq. (fT0]) , it is required 
for p s to be commutative to p, i.e. [p, p s ] = 0, which 
means that p s shares the same global symmetry with p. 
With this point one can determine g as Eq. (f2"Tj) . Fur- 
thermore since p s is diagonal under the collective basis 
{|n-f),77 = 0,l,--- ,N}, 

Tr[pp s ] = £ PnnP s nn < £ P "" + 2 (P "" )2 . (35) 

n n 

where p nn and p s nn denote the diagonal elements of p and 
p s respectively. Obviously the second equality occurs if 
and only if p nn = p s nn , which means that (J z ) has same 
value for both p and p s . And then a can be determined 
in Eq.pTj). 

As for the state 

|V) = -^(|1010> + 10101)), (36) 

our discussion above is also applicable. We should em- 
phasize our point clearly in this place that \ip) is not 
really translational invariance. Actually when one talks 
of the translational invariance about a systems, it means 

DHD^ = H, (37) 

in which D is the translation operator, and H is the 
Hamiltonian for this system. Hence that one speaks 
of the translational invariance for a state is meaning- 
less without specifying the Hamiltonian. Our discussion 
about Dicke model manifests clearly this point. So the 
crucial point is to find the Hamiltonian for which ip is 



the eigenvector. It seems that one can construct the fol- 
lowing Hamiltonian 

ff = XXo? +1 , ( 38 ) 

i 

for which ip is seemingly one of the degenerate ground 
states. If the translational invariance is required for this 
system, one must have the periodic boundary condition 
(7^r +1 = erf, where N is the total particle number. How- 
ever \ip) tells us that for one particle, its neighbored par- 
ticles always has opposite state to its state, which obvi- 
ously does not satisfy this periodic boundary condition. 
So we argue in this place that the translational invariance 
for I ip) is only occasional because of its special form. 
Instead \ip) is the true ground state for Hamiltonian 

ff=-E«+2> ( 39 ) 

i 

with the boundary condition o~ z N+2 — erf. It means that 
the particle always has the same state to its next neigh- 
bored particle. Thus it could explain naturally the reason 
that the maximal overlap with \ip) happens for the fully 
separable states 1 1010) or 1 0101) , which obviously satisfy 
this boundary condition and also are the ground states 
for this Hamiltonian. 

In another point, one can also find a state which seem- 
ingly satisfies the requirement of the "translational in- 
variance" defined by \ip), i.e., 

P = ^(|0101) (0101| + |1010)(1010|), (40) 

which obviously maximize the overlap with \ip). This 
features demonstrate again that \ip) is not truly transla- 
tional invariance since p' is the incoherent superposition 
of the two degenerate ground states for Eq. (|39| . 

Appendix B: Derivation of Eq. (f28)l 

Set 

Hq = uja^a; 

2A 

Hi = u J z + -=((J +a). (41) 
v A 

Under f3 = tt^f^ 1, the partition function can be ap- 
proximated as[43j 

z = Tr[e-^+^>] 

= Ti[e-^ Ho / 2 e- 0H ^ 2 e-^ H «/ 2 +0(/3 3 )} 

~ Tr[e-P Ho e~ pH ']. (42) 

With respect of 

N 

P s = £ N C n a n (l - a) N - n \n - f ){n - f |, (43) 
n=l 
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then 



A = -Ti-[pp s ] 



N 



= -Jr[J2NC n a n (l-a) N - n 

k=l 



(44) 



for which [p s , Ho] — is applied. The tricky for the 
tracing in the equation above is noting that | n — -y) 
denotes the state in which n particles are spin-up, and 
the others are spin-down. And then 



( n -%\ e -PB 0e -PB 1 \ n _N ) 



N 



'-" "/„ ^lJJ e - /3[i f' T? + ^ (at+Q)T 



N\ 



where L m (x) is the mth Laguerre polynomial, and the 
relation is used 



lm\(at + a) 2q \m) = -^ r {m\e r>{a +a) \m) 
drj z i 

d 2q -n 2 

' r e-L m (- V 2 )} 



r)=0 



(49) 



Apply the relation 

oo 

e- m ^L m (- V 2 ) 



exp[?7 2 



m=0 



and then 



1 - e-P" 



-I 



I, (50) 



oo olh 2 \2 
-/Wa/„ _ N\ ^ T^yl^ + ^(a* + a) 2 l fc 



2 I 



fc=0 



(2/c)! L 4 N 



A = I- 1 



z 1 - e-^ 



E 



{i 



Ah=0;fc2=0'**fcv=0 



{oo /q2A: 2 \ 2 



V -^-[^ + — fat + a) 2 l fc fl + £ 
A- (2fc)! [ 4 + A 1 + j M + 2fc + l 2 J 



ft go. 
2fc + 1 2 ' 

N —n 



n 



(2fci)! L 2fc + 

K=fci+fc 2 H hfcjv 

x E 



^ r f(l-2a)] 



,fc=0 



(45) 



9=0 

d 29 2- coth £_ 

e 2 2 



A 



Thus 



- f oo /j2fe 2 \2 



d?7 2 9 

With the relations 



(51) 



r/=0 



. fc=0 



[1 + 



w 



2fc + l 2 

Expand the item in the corbeil bracket 



d 29 22_ coth ^_ 
-e 2 2 



(l-2a)U ] (46) 



d?7 2 9 



(29- l)!!coth 9 ^ 



r;=0 



/ A />OC 

(2g- 1)!! = J-2 p A p / dxe^^x^, (52) 

V 7T J-oo 



E ^n^yyt 1 ■ 2fc + 1 2 



+ _J____(i _ 2a )] ] x and set A = 1/2, 



fci=0;fc 2 =0---fc 
7\=fci + fc 2 H \-k N 

E 

9=0 



a!(A-a)! 1 2 j W 1 + 



2 1 - e"/ 3 " 



fc 1 =0;fe 2 =0---fciv=0 



Define a'a|m) = m\m), and then 
1 



2 



N fl2fe, 



A = - V 

^ \ 11 (2fei)! L 2fc + l 2 

fe 1 =0:fe 2 =0---fcj V =0 \i=l v 1 
K=k 1 +k 2 A hfcjv 



(1 - 2a)] 



\1 = A ( 2 ^)! [ 2fc + l 2 

K=ki+ki-\ hfeiv 



(1 - 2a)] x 



E ^v(T) 2( "- 9) (^ C0t ^ 



9=0 



g!(A-a)! v 2 



N 



H 2 « »7 2 °° 

m=0 ^ = o 



Finally inverse the procedure from Eq. (|4"6f to Eq. (|4"7|) for 



the sum item and apply relations cosh a; 
(48)^inha; = ^—^ — , one then obtains the Eq. 



and 
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